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We study the interrelations between the classical (Frobenius-Perron) and the quantum (Husimi) 
propagator for phase-space (quasi-) probability densities in a Hamiltonian system displaying a mix 
of regular and chaotic behavior. We focus on common resonances of these operators which we 
determine by blurring phase-space resolution MM- We demonstrate that classical and quantum time 
evolution look alike if observed with a resolution much coarser than a Planck cell and explain how 
this similarity arises for the propagators as well as their spectra. The indistinguishability of blurred 
quantum and classical evolution implies that classical resonances can conveniently be determined 
from quantum mechanics and in turn become effective for decay rates of quantum correlations. 



I. INTRODUCTION 



Classical Hamiltonian motion can appear irreversible, as is most familiar for many-body systems but also observed 
for chaos of two degrees of freedom. Effective irreversibility can show up as relaxation of certain correlation functions, 
and the decay rates involved are resonances of the propagator for phase-space densities known as the Frobenius-Perron 
operator. Frobenius-Perron resonances have recently been shown to play a key role for the interrelations of classical, 
\^ | semiclassical and quantum dynamics, e.g. in a superanalytic approach to universal fluctuations in quantum (quasi-) 
£H ■ energy spectra ^-g].The purpose of this paper is to substantiate those interrelations. 

A quantum analog to the Frobenius-Perron operator is the propagator for Husimi functions. Husimi functions are 
positive and normalized distributions in phase space and turn into phase-space densities in the classical limit. Both 
classical densities and Husimi functions can be approximated within the same Hilbert space of phase-space functions. 
*y-\ • As a caveat we remark here that in the quantum case this Hilbert space of phase-space functions is not to be confused 
;> with the Hilbert space of quantum wave functions. 

Even though without coupling to external reservoirs both the classical and the quantum dynamics are reversible, 

both do appear effectively irreversible when phase-space structures cannot be fully resolved. Since any measurement of 

phase-space variables can be conducted with finite precision only, restricted phase-space resolution is rather ubiquitous 

in the macroworld; in particular for classical chaotic dynamics, where phase-space structures extend over an infinite 

y— ( . hierarchy of scales, but also in the quantum case, if a Planck cell is too small to be resolved. 

By adopting finite phase-space resolution much coarser than a Planck cell for a prototypical dynamics we show 
^h . that the quantum propagator becomes indistinguishable from the classical propagator. In particular, we find that 
resonances of the two propagators coincide, indicating that the same "relaxation processes" occur in classical and 
quantum dynamics. These processes and thus the resonances can be linked to the resolved structures in phase space 
Jupl. In particular, in as much as resonances are related to localized structures in phase space, our results suggest 
that resonances may cause some quantum eigenfunctions to assign exceptional weight to certain regions in phase or 
configuration space and possibly be one of the origins of scarring. This line of thought will be taken up in a separate 
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paper. 



If the classical phase space displays a mix of regular and chaotic structures, along with effective irreversible behavior 
caused by hyperbolicity, we also encounter recurrences due to elliptic islands in phase space; the latter manifest 
themselves as (almost) unimodular eigenvalues of the Frobenius-Perron and Husimi propagators. 

II. CLASSICAL DYNAMICS ON THE SPHERE 

As a prototypical Hamiltonian system with a mixed phase space we consider a periodically kicked angular momentum 
vector 

J = (j sin 9 cos ip, j sin 9 sin ip, j cos 9) (I) 

of conserved length j, also known as the kicked top. Such a system has one degree of freedom, and its phase space 
is the sphere, with the "azimuthal" angle ip as the coordinate and the cosine of the "polar" angle as the conjugate 



momentum. The dynamics is specified as a stroboscopic area-preserving map M on phase space. We choose the 
dynamics to consist of rotations R z (/3 z ),Ry(f3 y ) about the y— and z— axes by angles /3 y , Z and a "torsion", i.e. a 
nonlinear rotation T z (t) = R z (t cos 9) about the z— axis which changes <p by tcos9, 

M = T z ( T )R z (p z )Ry(p y ). (2) 

With P z and f3 y fixed, we can control the degree of chaoticity of the dynamics by varying r. While the dynamics is 
integrable for r = 0, with increasing r chaoticity sets in, until for r = 10 elliptic islands have become so small that 
they are difficult to detect. 

In the Liouville picture the time evolution of a phase-space density p is governed by Liouville's equation, 

d tP = Cp={H,p}, (3) 

where the Liouville operator C, the Poisson bracket with the Hamiltonian H, appears as the generator. For our 
rotations and torsion on the sphere we must separately take H — j3 y J y ,H — j3 z J z , and H — \J~l- Denoting the 
corresponding Liouvillians by C R , C Rz , Ct z we imagine Liouville's equation for each of them separately integrated 
over a unit time span. The product of the resulting three propagators yields the Frobenius-Perron operator V = 
expOCrJ exp(£flj exp(C Ry ). 

While a phase-space density is usually considered as L 1 — integrable, it can be assumed to belong to a Hilbert space 
of L 2 — functions as well. In this Hilbert space V is represented by an infinite dimensional unitary matrix. Accordingly 
its spectrum is unimodular and, depending on the character of the dynamics, may consist of both discrete eigenvalues 
and continuous parts. The complete basis to represent the matrix in may be chosen as ordered with respect to 
phase-space resolution, to eventually allow for a truncation of V to finite size in a systematic manner. 

On the unit sphere the spherical harmonics Yi m (6,tp) with I = 0,1,2,... and \m\ < I form a suitable basis. 
Employing that represenation resonances were identified in Ref . [HO] . We here briefly review those results that are of 
relevance for our present discussion. 

Phase-space resolution is characterized by the index I: By using all Yi m with < I < l ma x phase-space structure of 
area ex l/^ ax can be resolved. Our classical Liouvillians can be written in terms of the differential operators L y , L z 
well known from quantum mechanical contexts in the Y\ m representation, 

£r s = —^P Z L Z , L z — — i— — , 

dip 

( d d 

£ii y = -ip v L y , L y = -i I - cos ip — + cot 9 cos ip — 

C Tz = -ir cos 9L Z . (4) 

The Frobenius-Perron operator of the z— rotation becomes the matrix 

(exp(£ R j) = [d6sin6dpY l l(e > cp)exp(-(3 z -?-)Y l , m ,(0,cp) 

\ I Im.l'm' J \ Oip J 

= 5u<8 mm > exp(-im/5 z ) (5) 

which is diagonal in both indices. The Frobenius-Perron matrix of the y— rotation is blockdiagonal (diagonal in / but 
not in m) and consists of the Wigner d-matrices well known from quantum mechanics as 



(exp(£ fl J 



, „ = 5w& l mm ,{fly) . (6) 

Lrn.L' in' 



Finally, for the z— torsion the Frobenius-Perron matrix elements are finite sums over products of spherical Bessel 
functions ji(x) and Clebsch-Gordan coefficients, 



L-\-L 

(exp(£ T j) ,=<W(-irV(2*+l)(2Z' + l) £ (-if Ji"(mr)C l Zc u ;S . (7) 

V / LrnA'rn' — 

l" = \l-l'\ 

While both rotation matrices are diagonal in the indices 1,1', the torsion matrix T z is diagonal in m,m' . Therefore 
the elements of the Frobenius-Perron matrix for the composite dynamics are the products 

Vim,i> m > = ( exp(£ r J) (exp(£ijj) (exp(C Ry )) , (8) 

V / lm,l'm \ / l'm,l r m V / l'm,l r rn f 

and no infinite sum hinders their evaluation. 



III. QUANTUM EVOLUTION OF HUSIMI FUNCTIONS 

In the quantum mechanical description of the kicked top a wave vector \ip n ) is propagated by a Floquet operator 
F as IV'n+i) — F\il>n) over one period in between two kicks. Again F = T z (T)R z (/3 z )R y (f3 y ) consists of rotation and 
torsion operators that are built of components of an angular momentum vector J = (J x , J y , J z ) H as 

T z (t) = exp f-i - ■ + l Jz i . 
i? z (/3,) = exp(-i/3 z J z ) , 
R y ((3 v ) = exp(-ip y Jy) . (9) 



An obvious choice for the basis of the Hilbert space of wave functions are the (2j + 1) eigenvectors of J z , J z \jm) = 
m\jm) with fixed j and —j < m < j. While T z (t) and R z (f3 z ) are diagonal matrices in this representation, R y is 
again given by the Wigner d-matrix d 1 ((3 y ). Once again we would like to emphasize that the Hilbert space of quantum 
wave functions must not be confused with the Hilbert space of functions on the classical phase space employed in 
Sec. II. (The distinction between the two Hilbert spaces becomes especially important when one associates, for integer 
values of j, the vectors \jm) with the spherical harmonics Yj m and the components of J with differential operators.) 
As the quantum mechanical analogue to classical phase-space densities and the Frobenius-Perron operator we 
consider Husimi functions and their propagator. A Husimi function turns into a classical density in the classical limit. 
The Husimi function Q p of a density operator p is obtained as its diagonal matrix element with respect to a coherent 
state [0-0. Coherent states \j6p>) on the sphere assign to the observable J a direction characterized by the angles 
9 and ip, but contrary to the classical description that direction is only specified up to the minimum uncertainty 
permitted by the angular momentum commutation relations. The relative variance ((J 2 ) — (J) 2 )/i 2 = 1/j vanishes 
as the effective Planck constant 1/j goes to zero in the classical limit. Any coherent state \j9ip) can be obtained via 
a rotation R(9, ip) from the coherent state \j, m = j) as \j9<p) — R(9, p) \jj). The Husimi function thus reads 

Q p (9,p) = (j9v\p\j9p); (10) 

it can be expanded in terms of spherical harmonics, 

Qp = ^2 ^2 qimXlm , (11) 

1=0 m=— I 

with the index I limited to < I < 2 j [ Jl0[ . The finiteness of the latter expansion will be commented on presently. 

As a consequence, contrary to classical densities, Q p can be represented as a finite-dimensional vector with the 
dimension set by the effective (inverse) Planck constant j as (2j + l) 2 . Since no structures smaller than a Planck cell 
can be resolved, the Husimi propagator should be expected to be a finite matrix. 

The time evolution of Q p is governed by von Neumann's equation 

d t Q p (9, p) - QQ P (9, p) = -i{j0ip\ [H, p] \j9p) (12) 

with H = 2 ?_ ~ Jz,H = f3 z J z and H = (3 y J y for torsion and rotations, respectively. In the representation dll]) the 
generators Q can be written as differential operators very similar in structure to the classical Liouvillians (0) as (see 
Appendix A) 

/ 1 F) \ „ 

L z , 



Qt z = 


— it I cos 9 - 


1 d . 

-— sm 

2j + 1 89 


Gr s = 


-i-PzLz , 




Qr v = 


—i/3 y L y . 





(13) 

It is seen that the quantum Husimi generators Q and the classical Liouvillians C are identical for rotations, while for 
the z— torsion a quantum correction arises in Qt z , which is formally of order (2j + 1) . This correction is responsible 
for all the difference between classical and quantum mechanics in the sequel. We should point out here that the 



correction in question ensures the finite range of the Husimi matrix mentioned before, < I < 2j, — I < m < I. In 
contrast to the rotations, the torsion couples different angular momenta I. The action of the torsion generator on the 
spherical harmonics is calculated in Appendix B as 



.JUT (l + m+l)(l-m+l) v 

QTjlm = l^p^ I (2j - I) J (2| + 1)(2i + 3) Y l+1 , 



+<*+W&^H- (14) 

The generator couples spherical harmonics of neighbouring I, except for the cases I = and I = 2j, where the 
prefactors of Y/_i !m and Yz+i. m vanish, respectively. 

It follows from the foregoing remarks on the Husimi generator that the Husimi propagators exp(C^), exp(^j? ) for 
rotations are identical to their Frobenius-Perron correspondents. However, the torsion propagator exp(£/T z ) still has 
to be determined. Clearly, the (2j + l) 2 "Husimi functions" Q\j mi )(jm 2 \ = (j6 ( p\J m i)(j' m 2\j6 l p) are eigenfunctions 
of exp(Qr z ) with eigenvalues exp(— ir(mf — m\)/(2j + 1)). Note that the Q\j mi ){jm 2 \ represents the skew ket-bras 
\J rn i)(J rn 2 1 in the sense of the definition (flO|). In order to proceed to representing exp(C/T z ) in the basis {Yj m } we have 
constructed the linear transformation relating the spherical harmonics Yj TO to the "Husimi functions" Q\j mi )(jm 2 \ ( see 
Appendix C). Like its classical counterpart exp(£^) (see (M)) the Husimi propagator for z-torsion, exp(Cyj^), turns 
out diagonal in m but not in I. As a consequence, the Husimi matrix elements are simple products of torsion and 
rotation matrix elements, 

Uim,i> m > = ( exp(£ T J) (exp(g Rz )) (exp(^J) . (15) 

V / lm,l'm V / l'rn : l'rn V / l'rn,l'rn' 



IV. CLASSICAL TIME EVOLUTION AT FINITE PHASE-SPACE RESOLUTION 

As has been illustrated for the kicked top in two previous papers |y,|| , resonances and eigenvalues of the Frobenius- 
Perron operator as well as the associated phase-space structures can be found when time evolution is looked upon 
with finite resolution. However, since our aim is to compare classical and quantum dynamics at limited phase- 
space resolution, we briefly summarize the results. In representations with respect to the resolution-ordered basis 
functions Yi m with I = 0, 1, . . . , oo and integer \m\ < I, a limitation of phase-space resolution is achieved by a simple 
truncation, discarding all matrix elements of Vi m ,i'm' with I, I' > l max . The resulting N x N matrix V^ of dimension 
N = (lmax + I) 2 is a nonunitary approximation to the unitary V with a purely discrete spectrum bounded in modulus 
by unity. An investigation of the N— dependence of the eigenvalues of "p( N > reveals that with increasing resolution 
some eigenvalues persist in their positions, either close to or well inside the unit circle. Of these "frozen" eigenvalues 
the almost unimodular ones turn into unimodular eigenvalues of the unitary V as N — > oo; the subunimodular 
eigenvalues reflect resonances of V in a higher Riemannian sheet. The different nature of "frozen" unimodular and 
subunimodular eigenvalues becomes more evident through an investigation of their eigenfunctions. For unimodular 
eigenvalues the eigenfunctions are sharply localized on islands of regular motion surrounding elliptic periodic orbits, 
while for subunimodular eigenvalues eigenfunctions are localized around unstable manifolds of hyperbolic periodic 
orbits. The condition for freezing of an eigenvalue evidently is that at least the largest associated phase-space 
structures are resolved. For the kicked top with r = 10, f3 z — 1, and (3 y = 1 we propose to illustrate the freezing 
of eigenvalues. As no elliptic structures can be resolved even at the highest employed resolution, l max = 70, all 
frozen eigenvalues except for the one at unity (pertaining to the stationary constant eigenfunction loo) he well inside 
the unit circle. A gray-scale histogram of all eigenvalues with moduli greater than 1/4 of all matrices V^ with 
lmax = 20, 21 . . . , 70 in the complex plane is shown in figure |l|. Eigenvalues of smaller modulus have been rejected, 
since they have not settled yet and would spoil the histogram due to their large density near the origin. Dark areas 
in the histogram indicate large amplitudes and thus the positions of "frozen" eigenvalues. Several nommimodular 
"frozen" eigenvalues are clearly visible. Their precise positions at resolutions l max = 30,40,50,60, and 68 are given 
in table |. 



V. QUANTUM TIME EVOLUTION AT BLURRED PHASE-SPACE RESOLUTION 

The smallest scale in phase space resolvable in Husimi functions is set by Planck's constant. Recalling that for the 
kicked top the role of the inverse Planck constant is played by j, in the Hilbert space spanned by the Y/ m the Husimi 
propagator can be represented as a finite N x N matrix of dimension N = (2j + l) 2 . This is because with respect to 
the basis functions Yi m all matrix elements of U involving indices I > 2j vanish. The (2j + 1) 2 discrete eigenvalues of U 
are of unit modulus due to the unitarity of Schrodinger's time evolution which also implies conservation of probability 
in phase space. With the eigenvalue problem for the Floquet operator F\<j>k) = ex-p(—i<fik)\4>k) solved, the eigenphases 
of the full U are the quasienergy differences {<f>k — 4>i) and the associated eigenf unctions follow as (j9ip\(f)k)((j)i\j0<p)- 

We truncate the propagator U in the same manner as the Frobenius-Perron matrix, with the resolution parameter 
Imax taking a value l m ax < 2j. Now the dependence of the eigenvalues on the truncation parameter l max can be 
studied. We are particularly interested in semiclassical propagators, i.e. large values of j. Introducing the ratio 
K = l m ax/(2j), K = l represents complete resolution and k< 1 the limit with Planck cells far from resolved. 

Starting from k = 1 and slowly decreasing k, the eigenvalues of IA feel the truncation as a small perturbation at first 
that slightly shifts them inside the unit circle. With a further decrease of k, eigenvalues are spread over the entire unit 
disk. When k finally becomes small, the eigenvalues of the truncated propagator matrix coincide with the eigenvalues 
of the corresponding Frobenius-Perron matrix "p("> of the same dimension. This is illustrated in figure |2j, where the 
Husimi eigenvalues of the kicked top with r = 10, j3 z — j3 y — 1 for k — 1,0.5,0.16 and 0.08 are compared to the 
eigenvalues of the Frobenius-Perron matrix at l m ax = 32. For the sake of the comparison we kept l max = 32 fixed also 
for the quantum propagator and changed k by choosing different quantum numbers j = 16, 32, 100, 200. As the most 
striking aspect of the spectral coincidence at small k, resonances as well as eigenvalues of the full Frobenius-Perron 
operator can be obtained from quantum dynamics by looking for frozen eigenvalues of the truncated Husimi propagator 
at different values l ma x in the limit k <C 1. Furthermore, the coincidence of eigenvalues implies that also the associated 
eigenfunctions are in agreement. In particular, the eigenfunctions of the truncated Husimi propagator display the 
same strong scarring on unstable manifolds and elliptic islands in phase space as their classical counterparts. 

The most drastic difference between the two propagator matrices lies, of course, in the finiteness of the quantum 
propagator. Moreover, the matrix elements with I near the quantum cutoff 2j differ appreciably: It is illustrated 
in Appendix B that the quantum correction of the torsion generator in the basis of spherical harmonics is of the 
order l/(2j + 1) relative to the classical part. On the other hand, on large phase-space scales, i.e. for I <C 2j, the 
Frobenius-Perron and Husimi propagators become indistinguishable as k — > 0. An expansion of the Husimi propagator 
in powers of (2j + 1) _1 yields 

(exp(gj) =(exp(£)) +7r^(Y---) ■ ( 16 ) 

V / lm,l'm' V / lm,l'm' 2] + 1 V c — ' / lm,l'm' 

Assuming l max fixed, for j — ► oo the quantum correction of a matrix element decays as l/(2j + 1). Accordingly, in 
this case we expect the mean squared deviation between the matrix elements of the two propagators to asymptotically 
vanish as l/(2j + l) 2 . Figure || confirms that expectation. 

VI. CONCLUSION 

In conclusion, for a prototypical dynamical system we have formulated time evolution in terms of phase-space 
distributions, both classically and quantum mechanically. The comparison of the two respective propagators at 
limited phase-space resolution unveils the classical character of the quantum dynamics on large phase-space scales. 
Only when quantum coherences on the scale of Planck's constant are resolved, the peculiarities of quantum mechanics 
arise. Due to this fact Frobenius-Perron resonances can be identified from the quantum propagator. Moreover, 
quantum propagation is described by a finite matrix and therefore more easily treated than the classical counterpart 
which is an infinite matrix. Our results on classical and quantum signatures of resonances might become helpful in 
explaining phase-space localization (scars) of quantum eigenfunctions. 

We gratefully acknowledge support by the Sonderforschungsbereich 'Unordnung und grofie Fluktuationen' of the 
Deutsche Forschungsgemeinschaft and support of the Minerva foundation. 

VII. APPENDIX 



A. Husimi generators for rotation and torsion 

A normalized coherent state can be written as [|7| [l0[| 

j 
\j9p) = (1 + aa*)- j e aJ ~ \jj) = (1 + aa*)~ j ^ 



2j 



'-"j-rnW' 



(17) 



where a is the complex parameter a — tan | e lv and j± = J x ± i j y are the familiar ladder operators. In this manner 
a Husimi function (UQ) can be expressed as a function of a and its complex conjugate a* , 

Q p (a,a*) = (l + ao*)- !W (3J|e Q * i +pe oi -|Ji). (18) 

In order to derive the generators defined by von Neumann's equation (|12|) we start with the following two identities 



(jj\e a J + P J Z e aJ - \jj) = (jj\e a J + pe aJ - (J z - aJ-)\jj) 

j-a^){jj\e a * S +pe a J-\jj) 



Oj|e aV+ J zP e aJ -\jj) = (jj\(J z - a*J+)e Q * J + pe aJ ~\jj) 

d_ 
da 



= l7-T^-)0j|e aV+ pe aJ iJj> 



Thus, the generator for a rotation about the z-axis becomes 

^Gr.Qp = (1 + aa*r 23 {jj\e a ' J + \j z ,p] e aJ - \jj) 

Pz L J 



= (1 + act*) 



-2j 



j -a 



tf 







J -a- 



(jj\e a J+ pe aJ ~\jj). 



da* J \ da / 

It is easy to see that the differential operator commutes with the prefactor (1 + aa*)" 23 such that 

J z ^Q P ={a-^-a*^Q p . 

Replacing the complex variables a, a* by the angular coordinates 6, <p we finally get 



d__ . d 
da da* 



,d_ 
dip ' 



whereupon the generator becomes 



# 



■Gr z 



■ d f 

dip 



Some more effort is required for the torsion generator, 

i^l±lg Tz Q p = (1 + aa*)- 2 i(jj\e a ' J + 

T 

After using (O) twice one obtains 



J*,P 



,aj_ 



jj) ■ 



i^l±lg T _Q p = (l +aa *)- 2j 

T 



= (l+aa*)- 2j 
= (l+aaT 2j 



<9a" 



j - a 



da 



(jj\e a ' J+ pe aJ -\jj) 



j -a 

, d 
da 1 



dc 



-2ja~ 2ja*^ 

da da* 



(jj\e a J+ pe aJ -\jj) 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



The commutation of the left differential operator with the prefactor generates a further term, while the right one 
commutes as in (ElJ) , 

Again we replace the complex coordinates by the spherical coordinates, 

1 _2aa_ = l_aa 1=cos ^ 
1 + aa* 1 + aa* 

and the generator finally becomes 

^■=- t H-^rrlHi- (29) 



B. The torsion generator matrix 

We here calculate the action of the classical and quantum torsion generators on spherical harmonics. Spherical 
harmonics are separable into tp- and ^-dependent parts as 



W ^ ) = ^|^ r(OT( , )e .-, (30) 

where P™(cos6?) are the associated Legendre polynomials pj. At first we focus on the more complicated action of 
the generators Ct z and Gt z on the Pj m (cos0), i.e. for the classical propagator 

cos 9 Pp (cos 6) = zPr(z) , (31) 

and for the quantum correction (cf. (p9[)) 



- 2JT1W ™° p r(co S e) =-^ r± (z - (i - ^) A) p rw . (32) 

By using the following recursion formulae for the PJ n (z) H], 

(1 - ^ 2 ) Ap/"(z) = -teiH*) + (* + ™)^i(*) - (33) 

(21 + l)zPr(z) = (l-m + \)PrU + (l + m)P£ x {z) , (34) 

one arrives at 

for the classical generator and 

' ( z-{\- z 2 ) A) P™(z) = ___i___((/ + i )(/ _ m + l)P^(z) 1(1 + m )P{\(z)) (36) 



2 j + 1 

for the quantum correction term. The associated Legendre polynomials can now be replaced by the spherical harmonics 
(p0|). Taking into account L z Yi m = mYi m , we finally see the classical torsion generator to act on the spherical 
harmonics as 



L T _ l (l-m+l)(l + m + l) v (l-m)(l + m) v 

] ~~ Ylm ~ V (2Z + l)(2Z + 3) Y{l+1)m + V (2Z + l)(2;-l) r( '- 1)m ' (37) 



rm 



and the quantum correction as 



. Gt,-£t, „ l + l l (l-m + l)(l + m + l) „ I l (l-m)(l + m) „ 

1 rm lm ~ 2j + l\] (2Z + 3)(2/ + 1) Y{l+1)m + 2j + 1 \j (21 - l)(2l + 1) ( '~ 1)m ' W 

The expression for Qr^Yim is already given in (pf). The matrix elements of Lt z and Qt z in the basis of spherical 



harmonics can easily be read off from (pq) and (14), respectively. It is noteworthy that relative to the classical matrix 
element the quantum correction, to be read off from (pq), is of order l/(2j + 1). Therefore matrix elements of the 
quantum generator with I ~ 2j differ significantly from their classical counterparts; most importantly, the quantum 
correction manifests the finiteness of the expansion of Q at I = 2j, in accordance with the remarks after (pL4|) . 

C. Husimi torsion propagator 



In order to obtain the propagator matrix (|15]) the torsion propagator exp(^T 2 ) is to be written in the basis of 
spherical harmonics. Therefore we start from the eigenrepresentation of the propagator and apply the respective 
linear transformation. The basis of right-hand eigenfunctions consists of the (2j + l) 2 "Husimi functions" 

Q\jm 1 ){jm 2 \ = (J0<p\jmi)(jm2\j$<p) , (39) 

with eigenvalues exp(— ir(m 2 — m|)/(2j + 1)). 

The left-hand eigenfunctions lie in the function space dual to the space of Husimi functions. Functions from this 
space are the so-called P functions |j] -Q. For an operator p the P function P p (6, ip) is defined as the weight of the 
diagonal mixture with respect to coherent states, 

P=^T~J d6sin6 [ d<pP p (8,<p)\jdip)(j6<p\. (40) 

In particular, the functions P\jm 1 )(jm 2 \ are the left-hand eigenfunctions of exp(Crr, ) and form a basis of this dual space 
that is biorthonormal to the basis Q\j mi )(jm 2 \- Therefore, the eigenrepresentation of the torsion propagator is given 

by @ 

eXp(^? T J= Y^ IQ|jmi)<jm 2 |))exp I -IT ±- - 1 2 j ((P\j mi ){j m2 \\ , (41) 

mi,m 2 =—j 

wherein the double brackets refer to the scalar product of the classical Hilbert space. 

The coefficients qi m (mi, m 2 ) of the linear transformation connecting Husimi functions and spherical harmonics and 
the coefficients pim(iTii, fn^) of the inverse transformation can be obtained as Wp 10 



qi m (mi,m 2 ) = ((Yi m \Q\j mi )(jm 2 \)) = / d # dip sin 9 Y^Q\ jmi ) (j ma \(9, cp) , (42) 

Plm(mi,m 2 ) = ((*im|-P|jmi)<jm 2 |)) / d # d( /? siD -^ Y lm P \jmi)(jm 2 \ (#, <f) ■ (43) 

While the first integral yields the real coefficients (products of Clebsch-Gordan coefficients, incidentally 



qi m [mi,m 2 ) = d m , mi -*n 2 \ 4n{2l + V 



(I + m)\ y \J—miJ \j-Tri2s 

y? (-l) k+m (l + m + k)l(j + mi )[(j + k- m 2 )\ 

X 2^1 (1 _m - IcWhUm 4-lA\OA 4-h 4- m 4- -\M ' ^ ' 



(I - m - k)\kl(m + k)\(2j + k + m + 1)! 



the coefficients p/ m (mi, m-i) are proportional to the qi m (mi, m 2 ), 



^4n(2l + 1) 

Pirn\mi,m 2 ) = ——— qi m {mi,m 2 ). (45) 

(2j + l)qio,jj 



Actually, the above formulae are not particularly useful for numerical purposes as they involve an alternating sum 
of very large numbers. Easier to employ are recurrence formulae for the coefficients. One formula we used reads 



Ql(m-1 



i)(mi,m 2 ) = — -. r(w(j,mi)qi m (mi + I, mi) - w(i,-m 2 )g; m (mi,m 2 - 1) ) 

w{l, —to V / 



(46) 



wherein w(l,m) = y/Ul + 1) — m(m + 1). The recursion starts with 
qu(mi,m 2 ) 



A I a JWTVs ( 2 J) ! V(W / (j+m 1 )!(j-m 2 7 



(47) 



Coefficients with negative values of to can be obtained from ones with positive to through the symmetry relation 

qi { -m)(mi,m 2 ) = {-l) mi - m2 qi m (m u m 2 ) . (48) 

Finally, the propagator matrix in the basis of spherical harmonics is obtained as 



f exp(g T J 



= S n 



Im.l'm' 



E . / , 2mm i — m 2 \ 
qi m [mi,mi - m) pi> m (mi,mi -mjexp I -it — ) . (49) 

ii=max( — j-\-7n, — j) 
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("max — <J*J 


t"inax ^U 


^max — £**J 


*"max — OU 


''max — Oo 


0.8018 


0.8116 


0.8205 


0.8103 


0.8074 


0.7513 


0.7457 


0.7547 


0.7470 


0.7455 


-0.0041 
±i 0.7225 


-0.0063 
±i 0.7431 


-0.0123 
±i 0.7515 


-0.0079 
±i 0.7517 


-0.0027 
±i 0.7435 


-0.7391 


-0.7432 


-0.7475 


-0.7510 


-0.7427 


-0.6639 


-0.6766 


-0.6746 


-0.6869 


-0.6955 


0.6470 


0.6628 


0.6777 


0.6597 


0.6727 


-0.6218 


-0.6443 


-0.6188 


-0.6377 


-0.6336 


0.5462 


0.5968 


0.5847 


0.5889 


0.5786 


-0.5141 


-0.5406 


-0.5363 


-0.5481 


-0.5448 



TABLE I. Frozen nonunimodular eigenvalues of ~p( N > with r = 10, /3 X — /3 y = 1 at the resolutions l r . 



30, 40, 50, 60, 68. 




" 1 Re A, 1 

FIG. 1. For the top with r = 10, /3 Z — (3 y — 1 a histogram from the eigenvalues of p^' at resolutions lmax — 20, 21 ... 70 
in the complex plane has large amplitudes (black) at positions of frozen eigenvalues. Since the mean density of eigenvalues 
increases drastically near the origin, the disk with radius 1/4 is not included in the histogram. 
























© 














o£ / 


'«• <• 


























/ " 


^^^''o ^* 00 * 00 °° 


^ ^ 


\ 








/ ° 


e ° ^ + s+ ♦ ♦*£ * + % 


o°^ 


° » 


"o 






/ O<> 


o * / +0 ** + *H*VW 


*\ 


„v 


° G 






J °° 


+ °o fe o + 2 °/ + <>ef + 


,;' 


W 


■^ 























6 ^ * °^ ^*ijy»**** 






|«te«s<» 


















\ " * 






%<■ 


♦* 


/'= 


=32 



1 -1 




Re A. 



FIG. 2. Eigenvalues of the truncated Husimi (o) and Frobenius-Perron (+) operators (r = 10, (3 Z — 1, f3 y — 1) superimposed 
at lmax = 32 with the quantum number j taking the values 16 (full resolution), 32,100,200. While for j = 16 the Husimi 
spectrum is unimodular, for j — 200 classical and quantum eigenvalues are in good agreement. 
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C<1 
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<c 


0.01 
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0.001 
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0.0001 


1 « 


1e-05 




1e-06 




1e-07 




1e-08 
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FIG. 3. Mean squared deviation between Husimi and Frobenius-Perron matrix elements of the top (r = 10, /3 Z = l,/3 y = 1) 
with Imax = 4, 6, 8, 10, 12, 14, 16, 18, 20 as a function of j in a double-logarithmic plot. The numerical data suggest a power-law 
decay of the mean squared deviation as j — > oo. 
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